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Abstract 

Electronic conduction through single molecules is affected by the molecular electronic structure as well 
as by other information that is extremely difficult to assess, such as bonding geometry and chemical envi- 
ronment. The lack of an independent diagnostic technique has long hampered single-molecule conductance 
studies. We report simultaneous measurement of the conductance and the Raman spectra of nanoscale 
junctions used for single-molecule electronic experiments. Blinking and spectral diffusion in the Raman 
response of both /jara-mercaptoaniline and a fluorinated oligophenylyne ethynylene correlate in time with 
changes in the electronic conductance. Finite difference time domain calculations confirm that these corre- 
lations do not result from the conductance modifying the Raman enhancement. Therefore, these observa- 
tions strongly imply that multimodal sensing of individual molecules is possible in these mass-producible 
nanostructures. 
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The molecular- scale limits of electronic conduction are of fundamental scientific interest and 
relevant to future technologies. Our understanding of electronic conduction through single small 



molecules has grown dramatical 



mechanica 

HflQ 



break _ 
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unctionslllQ, 



ly in the last decade thanks to improved techniques ^ includin 



iflBflflyiiii 



3|] single-molecule transistors (SMTs)|l4, 
19], nanoparticle dimersl2oi. noise characterization [21], and thermopower 



measurement[l22']. A major complication in interpreting these experiments is the lack of local 
imaging or spectroscopic tools that can assess the environment and presence of the molecule of 
interest. 

Over the same period, single-molecule spectroscopies have progressed substantially. In partic- 
ular, surface-enhanced Raman spectroscopy (SERS) has been demonstrated with single-molecule 



sensitivity 1123 



25L 



26h in random aggregates of metal nanoparticles, though this level of de- 



tection is very challenging to prove conclusively. The electromagnetic component of SERS en- 
hancement results from the excitation of surface plasmons in the metal, leading to local fields at 
the molecule enhanced by a factor of g{uj) relative to the incident field. The Raman cross section 
is then enhanced by g(ujyg{uj')'^, where u and u' are the frequencies of the incident and Raman- 
scattered radiation, respectively. Electromagnetic SERS enhancements exceeding 10^^ are needed 
to approach single-molecule detection for many molecules without resonant Raman effects [|27ll . 
Additional "chemical" enhancement is also possible due to electronic interactions between the 
molecule and the metal. The vibrational modes observed in SERS reflect both the molecule and 
the local environment and conformation of the molecule on the metal surface. Using a sharp tip 



to provide a very 



molecule levelpS 



arge local field enhancement, Raman sensitivity has approached the single- 



290 . Very recently we found that nanoscale gaps between extended electrodes 



are very effective as extremely confined SERS hotspots and may be mass produced with relatively 
high yields. I30I1 

In this paper we use these nanoscale gap structures toperform simultaneous measurements o ' 
electronic transport and SERS. In many previous papers[4, 5,1^ 7, 8, 

\v\ [isi [19I it has been established that conductance in such structures is dominated by roughly a 
molecular volume. The conductance as a function of time is observed to correlate strongly with 
the SERS signal in 11% of the junctions measured. Conductance changes correlate with sud- 
den changes in the intensity of sets of Raman modes ("blinking") and with spectral diffusion of 
mode positions. The data suggest that both SERS and conductance changes are most likely due 
to changes in conformation and binding of an individual molecule. The combined data provide a 
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great deal of information about the effect of molecular orientation and environment on both con- 
duction and SERS, although a detailed understanding of this correlated information is indeed a 
very significant theoretical challenge. The most likely explanation for these results is that single- 
molecule multimodal sensing is possible. This combined measurement technique also opens the 
possibility of direct assessment of vibrational pumping and local heating in single-molecule elec- 
tronic transport. 

Our nanogap structures are fabricated on ~ 1 cm^ pieces of ra+-doped Si wafer with 200 nm 
of thermal silicon oxide. The structures are defined using electron-beam lithography and e-beam 
evaporation of 1 nm Ti and 15 nm of Au. The initial nanoconstriction structure consists of two 
large pads connected by a single constriction as shown in Figure 1 . The constriction is approx- 
imately 500 nm long and 100-180 nm wide. Liftoff is performed in acetone or chloroform and 
the samples are then cleaned of organic residue by 1 minute exposure to O2 plasma. For devices 
incorporating para-mercaptoaniline (pMA), samples are immediately placed in a 1 mM solution 
of pMA in ethanol. Samples are soaked in the pMA solution for 12 to 24 hours as pMA self- 
assembles on the Au surfaces, followed by an ethanol rinse to remove any excess pMA. We have 
also measured devices using a fluorinated oligophenylene ethynylene (F0PE)[j3iQ that possesses 
a distinctive Raman spectrum (see Supporting Information). FOPE was assembled via standard 
base deprotection[32] from a 0.25 mg/mL solution of the thioacetate form of the molecules in 1: 1 
ethanol-chloroform solvent, prepared under dry N2 gas. 

The nanoconstrictions are converted into nanogaps via electromigration [|33|l . a thoroughly 



studied[|34, 



3511 process dominated by momentum transfer from current-carrying electrons to 



mobile atoms in the metallic lattice. Electromigration has been used extensively to prepare 
electrodes for single-molecule conduction measurements, with typical yields of ~ 10-20% for 
tunneling gaps inferred to contain individual molecules (based on statistics on thousands of 
junctions) 141 15, , uM- Each gap is electromigrated with an 

automated procedure to form an atomic-scale constriction with a resistance of ~3 kTl, which is 
then allowed to break spontaneously pq] . While the atomic-scale details of each gap are different, 
gaps with measurable tunneling currents are formed routinely with high yield, and recent advances 
in it situ electron microscopy [|34i 13711 do, in principle, permit detailed structural examinations of 
the resulting electrodes. The migration and subsequent electrical measurements are performed in 
situ on the sample stage of the Raman measurement system, in air at room temperature. 

Electrical contact to the junction under test is made via micropositionable probes. One digital 
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lock-in amplifier (SRS SR830) is used to source 50-100 mV RMS at 200 Hz onto one pad, while 
the other pad is connected to a current-to-voltage converter (either SRS SR570 or Keithley 482). 
The AC current (oc dl/dV) and its second harmonic (oc cPl /dV"^) are measured with lock- in 
amplifiers, while the DC component of the current is sampled at 5.0 kHz. The unusually large AC 
bias (much larger than necessary to measure differential conductance alone) is required because 
of an unanticipated complication: the illuminated, molecule-decorated nanogaps can also exhibit 
significant DC photocurrents due to optical rectification (to be described in a separate publication). 
The large AC bias is needed so that the AC current is detectable without the DC current signal 
overloading the lock-in input stages. We find no evidence that the 100 mV RMS AC bias degrades 
the nanogap or the assembled molecules. 

Optical measurements are performed using a WITec CRM 200 scanning confocal Raman mi- 
croscope in reflection mode. Devices are illuminated by a 785 nm diode laser at normal incidence 
with a diffraction-limited spot. A lOOx ultra- long working distance objective leaves sufficient 
room for the micromanipulated electrical probes to be inserted between the objective and sample. 

As reported previously, the electromigrated nanogaps are outstanding substrates for SERS.||30(] 
Initially, spatial maps of the underlying Si of unmigrated nanogaps are obtained to facilitate cen- 
tering of the Raman microscope over the nanogap to within 100 nm. Spatial maps of the integrated 
molecular Raman signal after electromigration demonstrate the localization of the SERS hotspot 
(Fig. ID, E). Raman spectra are taken with 1 s or 2 s integration times while the microscope objec- 
tive is held fixed over the migrated junctions. Electrical measurements on unbroken constrictions 
under various illumination conditions demonstrate that heating of the electrodes due to the laser 
is not significant. The inferred change in the electrode temperature at ~ 0.5 mW laser power was 
less than 2 K (see Supporting Information). 

When the conductance of the junction drops below the conductance quantum. Go = 2e^//i, a 
tunneling gap is formed, and simultaneous conductance and Raman measurements are performed. 
In situ measurements of the optical response of nanogaps during migration are presented in Figure 
2. Even prior to complete nanogap formation, partially electromigrated junctions show SERS 
enhancement of the assembled molecules once the resistance exceeds about 1 k£L. The appearance 
of SERS indicates that the local interelectrode plasmon modes are now excitable. Crudely, this 
implies that over an optical cycle the junction acts more like a capacitor than a resistor; that is, 
the RC time constant of the nanogap is comparable to one optical period. For a 1 \Sl nanogap 
illuminated at 785 nm this implies an effective nanogap capacitance at optical frequencies on the 



4 



order of 10s of attofarads. 

The measured Raman signal strength scales logarithmically with the resistance of the gap un- 
til resistances exceed 1-10 MQ. At higher gap resistances, the Raman signal takes on a roughly 
constant value with sporadic changes (corresponding to SERS blinking events). This decoupling 
of electronic transport and SERS at low conductances is not surprising a priori, since tunneling 
conductances vary exponentially with gap size, while local plasmonic structure is less sensitive. 
Blinking events are often not correlated with further changes in junction resistance. This means 
that molecules are present in a region of strong Raman enhancement, while the molecular-scale 
tunneling volume dominating interelectrode conductance does not contain molecules that are de- 
tectably contributing to the Raman signal. We discuss this further below. 

However, in 17 of 120 junctions using ;?MA and 4 out of 70 junctions using FOPE, there are 
strong temporal correlations between the fluctuations in the nanogap conductance and changes 
in the SERS spectrum. This yield is quantitatively consistent with the yield of tunneling gaps 
containing single molecules inferred in single-molecule transistor measurements. Examples are 
shown in Figures 3 and 4 (with further examples in the Supporting Information). In Figure 3 A a 
simple positive correlation between Raman intensity and dl/dV is observed for all Raman modes. 
In Figure 3B another positive correlation between Raman intensity and differential conductance is 
observed in a different junction. In this case spectral diffusion of the Raman lines occurs but does 
not correlate significantly with the conduction. In both Figs. 3 A and 3B the amplitudes (count 
rates) of strong Raman modes have similar relative changes as the differential conductance. 

Figure 4 is an example of a more complicated relationship between the conductance and the 
SERS spectrum. While sudden changes in the Raman spectrum are correlated in time with changes 
in the measured conductance, some increases in Raman intensity correlate with increased conduc- 
tance, while others correlate with decreases in conductance. Additionally, changes in the mode 
structure of the Raman spectrum clearly correlate with changes in conductance. In region A the 
Raman spectrum and conductance are constant; when the Raman spectrum changes in region B a 
three-fold increase in conductance is observed, though overall Raman intensity changes only for 
certain modes. In region C the spectrum changes yet again, and while a weaker Raman spectrum 
remains, the conduction drops significantly. In region D the same mode structure seen in B returns 
and the conduction is similar to that seen in B as well. Regions D and E have positive correlations 
between Raman intensity and conduction with the lowest conduction observed between D and E 
where the Raman intensity is also lowest. At F a switch from positive to negative correlation be- 
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tween the Raman intensity and conduction occurs and carries over to regions G and H. At region 
I a small change in the mode structure is observed correlating with a switch to positive intensity- 
conduction correlations, continuing through regions J and K. In region L there is a final change 
in mode structure resulting in negative intensity-conduction correlations exemplified in the three 
conduction spikes that occur when the Raman spectrum disappears. 



It should be noted that our pMA spectra are typically dominated by the 62 symmetry modes 113 811 . 
as was seen previously OOll . This is not surprising, as it is well accepted that 62 symmetry modes 
experience additional "chemical" enhancements in comparison to ai symmetry modes. We often 
only observe the 1590 cm^^ ai symmetry peak and not the other expected oi mode at 1077 cm^^. 
Strong spectral diffusion in both molecules with shifts as large as ±20 cm^^ have also been ob- 
served, and are clear in Fig. 3b. This surely limits direct comparison to spectra reported elsewhere. 



However, the measured spectra are quite consistent with one another and are qualitative 



than those seen in "bare" junctions contaminated by physisorbed exogenous carbon BOh . 



y different 



One possible concern could be that changes in metal configuration at the junction are respon- 
sible for the fluctuations in tunneling conductance and SERS intensity. This scenario is unlikely 
for several reasons. First, tunneling conductances depend exponentially on gap geometry; while 
dl /dV could change by a factor of ten for a 0.1 nm change in gap separation, it is very unlikely 
that the electromagnetic enhancement would be as strongly affected. Second, it is not clear how 
metal rearrangement could explain the observed changes in Fig. 4; this would require that the gap 
itself alternately grow and shrink, with some changes in metal geometry giving large dl / dV fea- 
tures with small Raman effects and others vice versa. Finally, during the events in Figs. 3 and 4, 
the continuum emission at low wavenumbers observed previously due to inelastic light scattering 
from the metal electrodes [30] is constant in time (see Supporting Information). 

Another concern is that changes in tunneling conduction in one part of the junction may alter 
the plasmon mode structure and affect Raman emission from elsewhere in the junction. Such 
a scenario could lead to correlations like those in Figs. 3 and 4 even if conduction and Raman 
emission are not from the same molecule. Given that the interelectrode conductance affects Raman 
emission (Fig. 2), it is important to consider this possibility. We have performed finite-difference 
time domain (FDTD) simulations of the optical properties of such junctions to assess this issue, 
and the results {vide infra) effectively rule out this concern. While the finite grid size (1 nm) 
required for practical computation times restricts the quantitative accuracy of these calculations, 
the main results regarding spatial mode structure and wavelength dependence are robust, and the 
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calculated electric field enhancements are an underestimate^ 

Figure 5 shows a comparison of calculated extinction spectra that characterize the plasmonic 
mode structure of the gap structure shown, for various values of interelectrode conductance con- 
necting the source and drain at the indicated point. An analysis of the instantaneous charge dis- 
tribution associated with the plasmon resonances in Fig. 5 shows that negligible charge transfer 
occurs between the two electrodes for conductances smaller than Gq. The FDTD simulations 
show that the mode structure and enhancement are unaffected by conductances smaller than a few 
Gq. Details are presented in Supporting Information. These calculations confirm the interpreta- 
tion given above for Fig. 2: The plasmonic mode structure responsible for enhanced local fields 
in the nanogap is established once the interelectrode conductance falls well below Gq. Given 
these FDTD results, the only plausible explanation for such strong correlations in time between 
conduction and Raman emission is that both processes involve the same molecule or molecules. 

Conduction in nanogaps is known to be dominated by a tiny volume inferred in single molecule 
transport and break junction experiments to contain often only one molecule. Molecular move- 
ment, changes in bonding, or reorientation of the molecule in the gap results in different tunneling 
configurations and hence in conductance changes. The complex relationship between conduc- 
tance and Raman mode structure and intensity is also then natural, since chemical enhancement 
effects and the appearance of broken symmetry 62 modes should be strongly affected by changes 
in molecular configuration on the metal surface. It should also be noted that the measured junction 
conductances are consistent with the expected single molecular conductance range of 10^^ to 10~^ 
Go measured in similar molecules by break junction techniques [sj]. 

The detailed SERS mode structure combined with the conductance contains a wealth of in- 
formation about the bonding, orientation, and local environment of the molecule. With appro- 
priate electronic structure calculations and theoretical estimates of the Raman tensor for candi- 
date molecule/metal configurations, it should be possible to infer likely junction geometries and 
chemical structure corresponding to each type of Raman spectrum. Such calculations are very 
challenging even for the conductance distribution alone [|4o|. 

These conductance/Raman observations and accompanying calculations demonstrate that elec- 
tromigrated nanogaps between extended electrodes can achieve enhancements sufficient for single- 
molecule SERS sensitivity. Given that these structures can be fabricated in a scaleable manner in 
predefined locations with high yields OOll . this may allow significant advances in SERS-based sens- 
ing as well as multimodal sensing. With further improvements in technique (e.g., measurements 
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in vacuum as a function of temperature, interelectrode bias, and gate voltage in a SMT configu- 
ration), it will be possible to address open fundamental issues in SERS, including the nature of 
chemical enhancement, the mechanism of blinking, and the cause of the large spectral diffusion of 
Raman lines. Finally, comparisons of Stokes and anti-Stokes Raman peak intensities as a function 
of bias across the junction can reveal whether current flow pumps particular vibrational modes out 
of thermal equilibrium. This would enable new and detailed studies of nonequilibrium physics and 
chemistry at the single molecule scale. Comparisons between these results and transport/Raman 



studies on molecular ensembles [14 iL 



4211 should also provide valuable information about the effect 



of molecular environment on these processes. 
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FIG. 1: (A) Schematic of the electronic measurement. A 100 mV RMS AC signal is sourced by a lock-in 
into one pad. The AC current and its second harmonic are measured by lock-in amplifiers. The DC current 
is sampled with a current-to-voltage amplifier at 5 kHz. Raman spectra are synchronously captured with 
1-2 s integrations at an incident 785 nm wavelength laser intensity of about 0.5 mW. (B) Scanning electron 
image of Au constriction with nanogap. The constriction is 180 nm wide with a gap < 5 nm in size. (C) Map 
of the substrate Si 520 cm^^ peak (integrated from 480-560 cm^^) of the device from B. Red corresponds 
to the highest number of CCD counts and blue is the fewest counts. The Au pads which attenuate the Si 
signal are clearly visible. (D) Map of the pMA SERS signal from device in B from the ai symmetry mode 
at 1590 cm^^ (integrated from 1550-1650 cm^^), showing that the Raman signal is localized only to the 
nanogap region. (E) Map of integrated continuum signal (due to inelastic light scattering from the metal 
electrodes) from device in B (integrated from 50-300 cm^^). 
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FIG. 2: Blue curve (left scale): Resistance as a function of time for a nanogap as it is migrated. Migration 
is complete when the resistance reaches ~ 13 kO, Rr: (I/Gq)- Breaks in the curve occur where the gain 
of the current amplifier was being changed to maintain signal. Red curve (right scale): CCD counts per 
second in the 1590 cm^^ peak (integrated from 1550-1650 cm^^) as a function of time for the same device 
(synchronized with resistance plot). The intensity of the peak increases linearly with the log of the resistance 
until the resistance reaches around 10^ Q, at which point the intensity drops significantly and no longer 
shows correlations with the resistance. Stochastic intensity fluctuations ("blinking") are observed beyond 
this point. 
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FIG. 3: (A) Waterfall plot of Raman spectrum (2 s integrations) and positively correlated differential con- 
ductance measurements (dark blue = 50 counts; dark red =160 counts) for a pMA sample. All Raman 
modes that are visible exhibit this behavior as illustrated by the 1498 cm^^ and 1590 cm^^ modes. The 
1590 cm^^ mode has been shifted upward on the lower graph with the gray line indicating zero CCD counts. 
Vertical red lines indicate points of rapidly changing Raman intensity and conduction. Structure of pMA 
after self-assembly onto Au is shown below. (B) Waterfall plot of Raman spectrum (Is integrations) and 
positively correlated conductance measurement for a FOPE sample (dark blue = counts; dark red = 250 
counts). Strong spectral wondering is observed with no correlation to changes in conductance. Both visible 
modes at 1580 cm~^ and 1630 cm^^ are positively correlated. The slower response of the Raman spectrum 
compared to the conductance is due to the relatively long integration time. The 1630 cm^^ mode has been 
shifted upward on the lower graph for clarity. Structure of FOPE after self-assembly onto Au is shown 
below. 
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FIG. 4: Waterfall plot of Raman spectrum (Is integrations) and conduction measurements for a pVlK sam- 
ple. The device experiences periods of correlation (regions B,D,E) and anticorrelation (region L) between 
Raman intensity and conduction. Distinct changes in conduction are observed with every significant change 
in the Raman spectrum and are indicated by vertical red lines. The modes near 1394 cm^^ and 1522 cm~^ 
show similar intensity fluctuations except at region B and the end of region L. The color scale (dark blue = 
20 counts; dark red = 200 counts)has been set to make as many Raman modes visible as possible visible. 
This results in the saturation of the signal at region A which would otherwise resolve into well defined 
peaks. The 1522 cm^^ mode has been shifted upward on the lower graph for clarity. 
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FIG. 5: Extinction spectrum calculated using a 1 nm grid size for the structure partially shown in the left 
inset. The electrodes are modeled as Au, 15 nm thick, sitting on 50 nm thick Si02 dielectric, with an 
overall interelectrode gap of 8 nm. The upper and lower protrusions into that gap shown are modeled as 
hemispheres of radius 4 nm and 6 nm, respectively. The red square indicates the location of the modeled 
interelectrode conductance (a volume 2 nm on a side, meant to represent a molecule at the interelectrode 
gap). The right inset shows a map of |E|^, where E is the local electric field normalized by the magnitude of 
the incident field (roughly the Raman enhancement factor), for the mode near 825 nm. White corresponds 
to an enhancement of 10^. This field map is essentially unchanged until the junction conductance approachs 
10-^ S ~ Go. 
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Supporting Information: 
Simultaneous measurements of electronic conduction and Raman response in molecular 
junctions 

I. LASER HEATING 

To deteraiine the effects of laser heating on nanogaps, unbroken constrictions were illuminated 
with various laser powers while the resistance was measured as seen in Figure SI A. To determine 
the laser power we measured the number of CCD counts per second in the Si 520 cm~^ peak. 
At higher laser powers irreversible changes to the constriction occur as evident in the dropping 
resistance of the constriction when the laser power is held constant. We do not believe that such 
changes in the constriction will dramatically change our estimate of the laser heating. These laser 
powers are much higher than used in our conduction experiment where we expect about 50 CCD 
counts per second in the Si peak. From our measurements we determined the resistance of the 
constriction as a function of laser power as shown in Figure SIB using a linear least-squares fit. 
The final resistance/laser intensity relationship of = 1.35 x 10~^ Iiaser + 180 Q, where R is the 
resistance in Ohms and haser the laser intensity in Si CCD counts/second, is based on the average 
of three different samples. 

We then took different devices from the same sample chip and heated them in a probe station 
while measuring their resistance, as seen in Figure SIC. A linear fit was performed to determine 
the resistance as a function of temperature this time averaging over four samples. The final rela- 
tionship was determined to be = 0.1 IT + 85 Q, where R is the resistance in Ohms and T is the 
temperature in Kelvin. Equating our two relationships we find that the change in temperature per 
CCD count in the Si peak is 0.012 K. At the laser power used in our conduction experiments (50 
CCD counts/second) this lets us estimate a temperature change of 0.62 K due to illumination. 
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FIG. S 1 : (A) Resistance and laser intensity (in Si peak CCD counts/second) for an illuminated constriction 
as a function of time. It is clear that the resistance scales linearly with laser power. (B) Resistance as a 
function of laser intensity (again in Si peak CCD counts/second) for the same device as in A. The data is 
well represented by a linear fit of i? = 0.64 x 10~^ Iiaser + 170 Q where R is the resistance in Ohms and 
Iiaser is the lascr intensity in CCD counts/second. (C) Resistance as a function of temperature for a different 
constriction heated in a vacuum probe station. The data is fit to the line R = 0.085T + 92 i7 where R is the 
resistance in Ohms and T is the temperature in Kelvin. 
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II. ADDITIONAL EXAMPLES 



Figure S2 is a second example of the evolution of SERS response during the electromigration 
process. As in Fig. 1 of the main manuscript, Raman response becomes detectable at the nano- 
junction once the electromigrated two-terminal resistance exceeds 1-2 kTl. This implies that the 
localized interelectrode plasmon modes responsible for the enhanced local electric field are estab- 
lished before the junction is entirely broken. As discussed in the main text, a simple interpretation 
of this is that once the charge relaxation time across the nanoconstriction significantly exceeds the 
duration of an optical cycle, the coupling between the two sides of the junction is chiefly capacitive 
rather than conductive. 




Time [s] 

FIG. S2: Blue curve (left scale): Resistance as a function of time for a different nanogap as it is migrated. 
Breaks in the curve occur where the gain of the current amplifier was adjusted to maintain signal. Red curve 
(right scale): CCD counts per second in the 1590 cm^^ peak (integrated from 1550-1650 cm^^) of pMA 
as a function of time for the same device (synchronized with resistance plot). The intensity of the peak 
increases linearly with the log of the resistance until the resistance reaches around 10^ $7 where the Raman 
intensity then drops off significantly and no longer correlates with changes in the resistance. Stochastic 
intensity fluctuations ("blinking") are observed beyond this point. 
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Figures S3-S5 are additional examples of pMA Raman spectra that correlate in time with mea- 
sured junction conductance. Figure S6 shows similar data from one FOPE device and another 
pMA device. Sometimes the relationship between SERS intensity and G{t) is simple, as in Fig. 
S3A and S4A: increased Raman intensity correlates with increased conductance. However, in 
some devices decreases in SERS intensity, particularly of 62 modes, correlates with increased con- 
ductance, as in Figs. S3B and S4B. More complicated situations can also arise, as in Fig. S5, 
where conductance changes appear to correlate with both changes in Raman intensity and with the 
onset of spectral diffusion of some modes. These data contain much information implicit about the 
orientation of the molecule, its bonding to the metal atoms of the electrodes, and the local chemi- 
cal environment. A realistic theoretical treatment that encompasses these necessary ingredients is 
beyond the scope of this paper. 

III. CONTINUUM EMISSION AT LOW WAVENUMBERS 

As mentioned in the text, in these electromigrated junctions there is enhanced continuum Ra- 
man emission at low wavenumbers. This emission results from inelastic light scattering via the 
conduction electrons in the metal, and is present in junctions without molecules. Note that in de- 
vices with molecules, in extremely strong SERS blinks (brief periods of strong Raman emission), 
some increase in emission at low wavenumbers does occur, though this is much more modest than 
the percentage changes seen molecule-specific Raman modes. Figures S7 and S8 are reproduc- 
tions of the Raman data from Figures 3 and 4 of the main manuscript, but with the conductance 
as a function of time plotted in comparison to the integrated continuum emission (70 cm~^ - 
100 cm~^). During the major conductance changes and Raman blinking events, the continuum 
emission changes much less than either the SERS intensity or the conductance. This strongly sug- 
gests that the metal configuration at the junction is not changing in time; rather, the changes in 
SERS intensity and conductance are therefore due to changes in the molecular position, orienta- 
tion, and chemical environment. 
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FIG. S3: (A) Waterfall plot of Raman spectrum (2 s integrations) and positively correlated differential 
conductance measurements of a pMA sample (dark blue = 50 counts; dark red = 160 counts). Only the 
1525 cm^^ Raman mode exhibits a correlation with the conduction. The other modes as illustrated by 
the 1348 cm^^ mode do not show any changes expect at around 44 s where there is a temporary jump in 
the overall background. The 1525 cm^^ mode has been shifted upward on the lower graph with the gray 
line indicating zero CCD counts. Vertical red lines indicate points of rapidly change Raman intensity and 
conduction. (B) Waterfall plot of Raman spectrum (Is integrations) and negatively correlated conductance 
measurement for a different pMA sample (dark blue = counts; dark red = 250 counts). As illustrated by 
the modes at 1390 cm~^ and 1549 cm~^ conduction is only observed when the Raman intensity is almost 
undetectable. The 1549 cm~^ mode has been shifted upward on the lower graph for clarity. 
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FIG. S4: (A) Waterfall plot of Raman spectrum (1 s integrations) and positively correlated differential 
conductance measurements for a pMA sample (dark blue =10 counts; dark red =110 counts). For the most 
part the Raman intensity is very weak, however the intensity grows between 50 and 60 s culminating in a 
strong "blink" on and then off for the entire Raman spectrum. This is followed be an additional blink at 
around 70 s. Both the 1505 cm^^ and 1590 cm^^ modes show a correlation with the conduction during 
this period of blinking. The 1590 cm^^ mode has been shifted upward on the lower graph with the gray 
line indicating zero CCD counts. (B) Waterfall plot of Raman spectrum (Is integrations) and correlated 
conductance measurement for a different pMA sample (dark blue = 20 counts; dark red = 1000 counts). This 
sample experiences a more complicated correlations as seen in the at 1394 cm^^ and 1522 cm^^ modes. 
A spike in conductance is observed when the Raman spectrum completely disappears at around 18 s. The 
conduction and Raman intensity then fluctuate in a negatively correlated way. Finally at 30 s the conduction 
drops off while the mode at 1394 cm~^ continues to change. The 1522 cm^^ mode has been shifted upward 
on the lower graph for clarity. 
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FIG. S5: Waterfall plot of Raman spectrum (1 s integrations) and conduction measurements for another 
pMA sample. This device shows that the conduction may be correlated with spectral wandering rather 
than intensity changes in the Raman spectrum as is evident in the 2128 cm^^ mode. We observe spectral 
wandering changes between 25 and 55 s that appear to correlate with alterations in the conduction. At 160 s 
a drop in the overall background is observed that artificially lowers the intensity of all modes but has no 
effect on conduction. At 213 s the mode at 2128 cm~^ along with others vanishes and the conduction also 
drops significantly and returns along with the Raman modes at 220 s. Finally a strong shift in wavenumber 
for the peak at 2128 cm^^ is observed at 265 s corresponding with a drop in conduction at the same time. 
The 2128 cm~^ mode has been shifted upward on the lower graph for clarity. 
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FIG. S6: (A) Waterfall plot of Raman spectrum (1 s integrations) and positively correlated differential con- 
ductance measurements for FOPE molecule. The Raman intensity is plotted on a log scale for clarity(dark 
blue = 4 counts; dark red = 410 counts). The Raman intensity for the 350 cm~^ and 1077 cm^^ is observed 
to steadily peaking at 62 s for both modes and again at 80 s for the 1077 cm^^ mode. The conductance 
also peaks at 62 s but also has some obvious differences in hne shape from the Raman intensity. However, 
the log of the conduction shows good correlation with the 350 cm~^ peak. The 1077 cm^^ mode has been 
shifted upward on the lower graph with the gray line indicating zero CCD counts. The slower response 
of the Raman spectrum compared to the conductance is due to the relatively slow integration time. (B) 
Waterfall plot of Raman spectrum (Is integrations) and positively correlated conductance measurement 
for a pMA device (dark blue = 60 of fewer counts; dark red = 200 counts). Strong spectral wondering is 
observed with no correlation to changes in conductance. Both visible modes at 1160 cm^^ and 1540 cm^^ 
are positively correlated, no al symmetry modes are visible. The slower response of the Raman spectrum 
compared to the conductance is due to the relatively slow integration time. The 1540 cm^^ mode has been 
shifted upward on the lower graph for clarity. 
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FIG. S7: Reproduction of figure 3 from original paper. Here the low wavenumber continuum intensity 
(integrated from 70-100 cm~^) is plotted as a function of time. In both A and B the continuum is observed 
to be relatively constant with fluctuations much smaller than those seen in the Raman modes. 
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FIG. S8: Reproduction of figure 4 from original paper. Here the low wavenumber continuum intensity 
(integrated from 70-100 cm~^) is plotted as a function of time. The continuum shows little correlation with 
the overall conduction especially when compared with the strong correlation exhibited by the active Raman 
modes. The observed jumps in continuum typically correspond with times when the background of the 
entire spectrum is elevated relative to very high wave numbers (> 3000 cm^^). This could be the result of 
continuum generation from the molecule. 
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IV. FOPE MOLECULE 



Figure S9 shows the Raman spectrum for a bulk crystal of the FOPE molecule, first reported 
by Hamadani et a/. [SI] The mode at approximately 2210 cm~^ is the symmetric stretch of the 
carbon-carbon triple bond. The full name of FOPE is thioacetic acid 4-pentafluorophenylethynyl- 
phenyl ester, and its synthesis and characterization are reported in the supplemental material to the 
above paper. 
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FIG. S9: Raman spectrum of bulk FOPE crystal taken in Renishaw Raman microscope at 785 nm. Curvature 
and offset of spectrum is due to fluorescence. 
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V. FDTD CALCULATIONS 



The optical properties of the bowtie structure were calculated using the Finite-Difference Time- 
Domain method (FDTD) using a Drude dielectric function with parameters fitted to the experi- 
mental data for gold. This fit provides an accurate description of the optical properties of gold for 
wavelengths larger than 500 nm.[S2] These calculations do not account for reduced carrier mean 
free path due to surface scattering in the metal film, but such an effect is unlikely to change the 
results significantly. 

The bowtie is modeled as a two finite triangular structures as partially illustrated in Fig. 5, left 
inset, of the manuscript. Our computational method requires the nanostructures to be modeled to 
be of finite extent. The plasmon modes of a finite system are standing modes with frequencies 
determined by the size of the sample and the number of nodes of the surface charge distribution 
associated with the plasmon. For an extended system such as the bowties manufactured in this 
study, the plasmon resonances can be characterized as traveling surface waves with a continuous 
distribution of wavevectors. 

Previous calculations of a series of calculations of bowties with increasing length reveals that 
the optical spectrum is characterized by increasingly densely spaced plasmon resonances in the 
wavelength regime 500-1000 nm and a low energy finite-size induced split-off state involving 
plasmons localized on the outer surfaces of the bowtie. For a large bowtie, we expect the plasmon 
resonances in the 500-1000 nm wavelength interval to form a continuous band.[S3,S4] 

To investigate the effects of a conductance shunting the nanoscale gap, FDTD calculations were 
performed for a bowtie with two semi-spherical protrusions in the junction as shown in the close- 
up in the inset of Fig. 5 of the main text. The electrodes are modeled as regular trapezoids of a 
height of 50 nm, and a 1 nm grid size was used. The conductance was modeled as a cubical volume 
2 nm on a side located between the 6 nm radius asperity and the facing electrode, where the local 
field enhancement is maximized for modes relevant to the wavelengths used in the experiment. 
The conductivity of the material was set to be frequency independent over the wavelength range 
of interest (as expected for tunneling), and chosen such that the conductance of that interelectrode 
link was the desired value. 

Figure SlO(top) reproduces the main portion of Fig. 5 from the main text. The three most 
prominent features in the calculated extinction spectrum are labeled. Peak "a" corresponds to the 
mode shown in Fig. 5 of the main text, believed to be most relevant for the experiments at hand. 
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Figure SlO(bottom) shows the evolution of the electric field enhancement factor (calculated for the 
mid-point of the conducting volume standing in for the molecule) as a function of interelectrode 
conductance. The field enhancement and mode shapes due to local features in the junction are 
essentially unaffected by the interelectrode conductance until that conductance exceeds the order 
of Go = 2e^ /h = 7.74 x 10~^ S. For conductances significantly larger than Gq, charges can 
flow between the two electrodes and a new low energy plasmon resonance appears at wavelengths 
that depend on the conductance of the shunt. This is completely consistent with the observations 
in Fig. 2 of the main paper and Fig. SI of this supporting material. This demonstrates that the 
conductance changes in these experiments are not able to modify Raman emission from other 
molecules not in the tunneling region. 
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FIG. SIO: Top: Extinction spectra via FDTD calculations as explained in the main text. Bottom: Field 
enhancement (relative to incident field) calculated at the midpoint of the interelectrode conductance region 
for the three modes labeled in the upper graph, as a function of interelectrode conductance. Very little 
change in field enhancement or distribution is seen until interelectrode conductance exceeds Gq. 
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